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ABSTRACT 

Context. Observed X-ray spectra of some isolated magnetized neutron stars display absorption features, sometimes interpreted as ion 
cyclotron lines. Modeling the observed spectra is necessary to check this hypothesis and to evaluate neutron star parameters. 
Aims. We develop a computer code for modeling magnetized neutron star atmospheres in a wide range of magnetic fields (10'- - 10'^ 
G) and effective temperatures (3 x 10^ - 10' K). Using this code, we study the possibilities to explain the soft X-ray spectra of isolated 
neutron stars by different atmosphere models. 

Methods. The atmosphere is assumed to consist either of fully ionized electron-ion plasmas or of partially ionized hydrogen. Vacuum 
resonance and partial mode conversion are taken into account. Any inclination of the magnetic field relative to the stellar surface 
is allowed. We use modern opacities of fully or partially ionized plasmas in strong magnetic fields and solve the coupled radiative 
transfer equations for the normal electromagnetic modes in the plasma. 

Results. Spectra of outgoing radiation are calculated for various atmosphere models: fully ionized semi-infinite atmosphere, thin 
atmosphere, partially ionized hydrogen atmosphere, or novel "sandwich" atmosphere (thin atmosphere with a hydrogen layer above 
a helium layer). Possibilities of applications of these results are discussed. In particular, the outgoing spectrum using the "sandwich" 
model is constructed. Thin partially ionized hydrogen atmospheres with vacuum polarization are shown to be able to improve the fit 
to the observed spectrum of the nearby isolated neutron star RBS 1223 (RX 11308.8+2127). 

Key words, radiative transfer - methods: numerical - stars: neutron - stars: atmospheres - X-rays: stars - stars: individual: RX 
J1308.8-H2127 



1. Introduction 

In the last two decades, several new classes of neutron stars 
(NSs) have been discovered by X-ray observatories. They in- 
clude X-ray di m isolated NS s (XDINSs), or Magnificent Seven 
] (see review by iHaberll 2007h. central com pact objects (CCOs) 
■ in supernova remnants (iPavlov et al.l2002l 2004), anomalous X- 
' ray pulsars and soft-gamma repeaters (AXPs and SGRs; see re- 
views by Kaspi 2007; Mereghetti et al. 2007; Mereghetti 2008]). 
The NSs in the last two classes have superstrong magnetic fields 
(B ^ 10'"* G) and are commonly named magnetars. The XDINSs 
can have B ~ a few x 10'^ G, as evaluated from period changes 
and from absorption features in the observed spec tra, if they are 
interpreted as ion cyclotron li nes (see reviews bv lHaberlll2007t 
Ivan Kerkwiik & Kaplanll2007l) . 

These NSs are relatively young with ages < 10^ yr and suf- 
ficiently hot (reff ~ 10* - 10' K) to be observed as soft X-ray 
sources. The thermal spectra of these objects can be described by 
blackbody spectra with (color) temperatu res from 40 to 700 eV 
(see, for example. iMereghetti et 111120021) . Some of the XDINSs 
and CCOs in supernova remnants have one or more absorption 
features in their X-ray spectrum at the energies 0.2 - 0.8 keV 
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dHaber]l l2007^. The central energies of these features appear 
to be harm onically spaced (Sanwal et al. 2002; Schwope et aT] 
I2OO7I: lyarTKerkwiik & KaplanI 120071: lHaberlll2007i) . The optical 
co unterparts of sorn e XDINSs are also known (see review by 
Mig nani et al.l2007h . and their optical/ultraviolet fluxes are a few 
times larger than the blackbody extrapolation of the X-ray spec- 
tra (Burwitz et al. 2001, 2003; Kaplan et al. 2003; Motch et al] 
l2()03h . 

The XDINs are nearby objects, and parallaxes of some of 
them can be measured (Kaplan etal. 2002a). Therefore, they 
give a good possibility to measure the NS radii, yielding useful 
information on the equation of state (EOS) for th e NS inner core 
dTrumper et al.ll2004t iLattimer & PrakashlllOOTi) . which is one 
of the most important problems in the NS physics. For example, 
the EOS is necessary for computations of templates of gravi- 
tatio nal wave signals w hich arise during neutron stars merging 
(e.g. lBaiottietal.ll2008l) . 

For a sufficiently accurate evaluation of NS radii, a 
good model of the NS surface radiation for the observed 
X-ray spectra fitting is necessary. The isolated NS surface 
layers can either be condensed or have a plasma enve- 
lope dLai & Salpeteij 119971; iLail I2OOII) . In the latter case, the 
outer envelope layer forms an NS atmosphere. The structure 
and emergent spectrum of this atmosphere can be computed 
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by u sing stellar model atmosphere methods (e.g. iMihalasI 
Il978l) . Such modeling has been performed by many scien- 
tific groups for NS model atmospheres without magnetic field 
( Romani 198 7; Zavlin^et al. 1996; Raiasocal & j^omani 199; 
IWerner & Deetienll2000l; lOansicke et al.H2002l; iPons et al.ll20a 
and by several groups for models with strong ( B <: lO'^ G) mag- 
netic fields jS hibanov et al. 1992; Raiaaopal etal.lll997t lOzell 
1200 U iHo & Laill200k i2003l l2004t |vanAdelsberg & Lai"2006|). 
These model spectra were u sed to fit the observed isolated NS 
X-ray spectra (see review bv lZavhnll2009l) . 

Modeling of magnetized NS star atmospheres is based on 
the theory of electromagnetic wave propagation in a magne- 
tized plasma in t wo normal modes, extraordinary (X) and or- 
dinary (O) ones (lGinzburglll970l; lMeszaroslll992l) . an d on the 
methods of opacity calculations for these two modes (IVentural 
119791; iKaminker et aTl Il982l [19831) . Methods of fully ionized 
model atmo spheres modeling are well developed (see, e.g., 
IZavlinI l2009l for references). P artially ionized hydrogen atmo- 
spher es have been modeled dPotekhin et al.l 120041; iHo & Lail 
2004 IHo et al. 20081). using t he opa city and EOS calculations by 
Potekhin & Chabrierl 12003', 2004), which accurately take into 



account the motional Stark effect and plasma nonideality ef- 
fects in quantizing magnetic fields. Mid-Z element atmospheres 
for str ongly magnetized NSs have been modeled bv lMori 'feHol 
(I2007h . who treat the motional St ark effect using a perturbation 
theory dPavlov & Meszaroslll993h valid at relatively low T. 

For magnetar atmospheres, polari zation of the vacu u m can 
be significant, which was studied by iPavlov & GnedinI (1 19841) 
and recently bv lLai & Hol ( l2002ll2003l!) . Model atmospheres with 
partial mode conv ersion due t o the vacu um polarization have 
been c omputed by iHo &Lail (I2003h and Ivan Adelsberg & Lail 
(Hop). 



If the temperature is sufficiently low or the magnetic 
field is sufficiently strong, th e thick atmosphere can be re- 
placed by a condensed surface dLai & Salpeted[T997t lLaill2001l; 
[Medin&Lai 200^ as a result of the plasma phase tran - 
sition (cf. Potek hin et all 119991: iPotekhin & Chabried |2004|) . 
Emission and absorption properties of such surfa ces i n stron g 
magnetic fields have been studied by [X urolla e t al.l (|2004|) ; 
Ivan Adelsberg et al.1 (l2005l) ; lPerez-Azonn et al.i 62005^ 

In recent years, evidence has appeared that some of the 
XDINSs may have a "thin" atmosphere above the condensed 
surface. Such atmosphere could be optically thick to low -energy 
photo ns and optically thin to high-energy photons. Mo tch et al.l 
(l2003l) fitted the spectrum of RX J0720.4-31 25 using a non- 
magnetic thin hydrogen atmosphere model, and lHo et al.l d2007l) 
fitted the spectrum of RX J 1856.4-3754 using a magnetic, par- 
tially ionized, thin hydrogen atmosphere model. 

In this paper we present a new computer code for model at- 
mospheres of magnetized NSs and some new results, which were 
obtained by using this code. The code is applicable to modeling 
of fully ionized NS atmospheres and partially ionized hydrogen 
atmospheres accounting for the partial mode conversion due to 
the vacuum polarization. The magnetic field can be normal or 
inchned to NS surface. We use this code to study various thin 
atmospheres above a condensed NS surface. In particular, we 
consider a novel "sandwich" atmosphere model, where a thin at- 
mosphere is composed of a hydrogen layer above a helium layer 
We also discuss the possibility to apply our thin partially ionized 
hydrogen model atmospheres with vacuum polarization to fit the 
spectrum of RBS 1223. Our numerical method is outlined in §|2j 
Results are presented and discussed in §|3] Conclusions are sum- 
marized in §|4] 
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Fig. 1. Geometry of the radiation transfer. 



2. Method of atmosphere structure calculations 

We compute model atmospheres of hot, magnetized NSs subject 
to the constraints of hydrostatic and radiative equilibrium assum- 
ing planar geometry and homogeneous magnetic field. There are 
two versions of the code. In the first one, we consider the mag- 
netic field B perpendicular to the surface. In this case the angle 
a between B and a radiation wave vector k is equal to the angle 
6 between k and the normal n to the surface (see Fig. [T]i. It is the 
simplest case, because opacities depend on a, and the geome- 
try of radiation propagation depends on 6. In the second version, 
the angle 9b between B and n is arbitrary, and calculations are 
more expensive. In this case the opacities depend not only on 
the polar angle 9, but also on the azimuthal angle (f between the 
projections of B and n onto the stellar surface, therefore it is 
necessary to solve radiation transfer equations for a significantly 
larger number of directions. 

The model atmosphere structure for a NS with effective tem- 
perature Teff, surface gravity g, magnetic field B, and given 
chemical composition is described by the following set of equa- 
tions: 



1 . The hydrostatic equilibrium equation 

AP„ 

-j^ = g- gmd, 
am 

where 

GMns 



4s VI -Rs/R 



NS 



^lad 



-j^ 2 poo p27T p+l 



(1) 



(2) 



(3) 



/■=! 



allows for the radiation pressure, and - IGM^^jc^ is the 
Schwarzschild radius of the NS. Here ji - cosO, I\,(p.,ip) is 
the specific intensity in mode /, is the gas pressure, and 
the column density m is determined as 



dm = —pdz. 



(4) 



The variable p denotes the gas density and z is the vertical 
distance. Of course, radiation pressure is unimportant in the 
models presented below, but it can be more significant at 
higher effective temperatures. 
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2. The radiation transfer equations for the two modes 
dp 

where 



" K + cr\,2 



(5) 



(6) 



\ \ ^ r^" r"^' 

ITTi tYj I ^"f' I crv'(/^.^;/^'.^')/v(y«',i^')dM' 

2n k\, + a-'yj-f Jo J-i 

is the source function. By is the blackbody (Planck) intensity, 
and the optical depth tJ, is defined as 



dT; = (fe; + Odm. 



(7) 



Here, the true absorption and electron scattering opacities k[, 
and (Ty depend on jj. and ip. Specific intensity in given direc- 
tion and mode can be scattered in some other direction and 
in both modes 



1 \' 1 C^^ 



(8) 



./=i 

3. The energy balance equation 



1 ^ /^oo r*ln 



(9) 



with emissivity 



2 



(10) 



Equations ([U - (fTOl i must be completed by the EOS and the 
charge and particle conservation laws. In the code two different 
cases of these laws are considered. 

In the first (simplest) case, a fully ionized atmosphere are 
calculated. Therefore, the EOS is the ideal gas law 



(11) 



where ntot is the number density of all particles. Opacities 
are calculated in the s ame way as in the paper by 
Ivan Adelsberg & Lail ( l2006l) (see references therein for the back- 
ground theory and more sophisticated approaches). The vac- 
uum polarization effect is taken into account following the same 
work. Examples of opacities as functions of photon energy and 
angle 6 for a magnetized electron-proton plasma are shown in 
Figs.|2]and[3l 

The second considered case is a partially ionized hydrogen 
atmosphere. In this case the EOS and the corresponding opac- 
ities are take n from tables calculated by Potekhin & Chabrie^ 
(I2003L |2004|) . The nor mal mode polarization vectors are taken 
from the calculations bv lPotekhin et al.l (|2004|) . The vacuum po- 
larization effect is also included. 

For solving the above equations and computing the model 
atmosphere, we used a version of the computer code ATLAS 
(iKuruczl 1 19701 Il993h . modified to deal with strong magnetic 
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Fig. 2. Dependence of the free-free and electron scattering opac- 
ities in two modes on the photon energy at different angles 
in a fully ionized hydrogen plasma. The proton cyclotron and 
vacuum resonances are also shown. The plasma temperature is 
7x10*' K, the plasma density is 30 g cm"-', the magnetic field 
strength is 10"* G. 



fields. A nonmagnetic version of this modified code was pre- 
viously used t o mod e l atmospheres of super- soft X-ray sources 
(ISwartz et al.l 1200 2': Ibragimov et al. 2003), atmosphe res of 
non-magnetized NSs (Suleimanov & Werner 2007; Rau ch et al.l 
2008), and at mospheres of spreading layers on the surface of ac- 
creting NSs ( Suleimanov & Poutanenl2006l) . 

The scheme of calculations is as follows. First of all, the 
input parameters of the model atmosphere are defined: T^s, g, 
B and the chemical composition. Then a starting model using 
a grey temperature distribution is calculated. The calculations 
are performed with a set of 90 depth points mj distributed log- 
arithmically in equal steps from m\ » 10"^ g cm"^ to m^^^ » 
10^ g cm"^ in the case of a semi-infinite atmosphere. It is also 
possible to calculate thin atmospheres with arbitrary values of 
Wmax- In this case the temperature at the inner boundary is con- 
sidered as the temperature of a condensed NS surface. 

In the starting model, all number densities and opacities at all 
depth points and all photon energies are calculated. We use 200 
logarithmically equidistant energy points in our computations in 
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Fig. 3. Dependence of the free-free and electron scattering opac- 
ities in two modes on the angle between the magnetic field lines 
and the direction of photon propagation for two photon energies: 
1 keV and at the proton cyclotron energy. Plasma parameters and 
magnetic field strength are the same as in Fig.|2] 
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Fig. 4. Dependence of the true opacities in two modes on the 
photon energy for fully (dashed curves) and partially ionized 
(solid curves) hydrogen. The vacuum polarization effect is taken 
into account. The plasma temperature is 10^ K, the plasma den- 
sity is 1 g cm"^, the magnetic field strength is 4 x lO''' G. The 
angle between photon propagation direction and the magnetic 
field is 60°. 



the range 0.001 -20 keV with 9 additional points near each ion 
cyclotron resonance 



£ci = 0.635 keV 



B 



1014 G A' 



(12) 



where Z is the ionic charge and A the atomic weight in the atomic 
units. If the vacuum resonance is taken into consideration, then 
another photon energy g rid is used, whic h is constructed using 
the "equal grid" method jHo & Laill2003l) . In this method every 
point in the depth grid mj corresponds to the point in the energy 
grid defined by the equation 



Ei = Ev(p(mj)) 



(13) 



(if this energy point is in the considered energy range 0.001 -20 
keV). Here E y(p) is the energy of the vac uum resonance at given 
B and p (see Ivan Adelsberg & Laill2006l) . The opacity averaged 
in the energy interval (£,-1, £,+1) is used to avoid opacity over- 
estimation in frequency integrals like (|9]l. This energy grid is re- 
calculated after every iteration. Note that the vacuum resonance 
energy Eyip) in a partially ionized plasma is shifted relative to 
its value in a fully ionized plasma (see Fig.|4]i. 

The radiation transfer equation (|5]l is solved on a set of 
40 polar angles 6 and 6 azimuthal angles ip (in the case of 
inclined magnetic fiel d) by the short characteristic method 
(l01son&Kunasdll987h . 

We use the conventional condition (no external radiation) at 
the outer boundary 



I[,(fi <Q,m- mi) = 0. 



(14) 



The diffusion approximation is used as the inner boundary con- 
dition 



fyifl > 0,m = mmax) = 



(15) 



The code allows one to take into account the part ial mode 
conversion according to Ivan Adelsberg & Lail ( 1200 6*). At the 
vacuum resonance, the intensity in the extraordinary mode par- 
tially converts with the probability 1 - Pjump to the intensity in 
the ordinary one and vice versa 



p. rl.2 /I _ p. w2,l 
'v jump ~ \^ ^ jump/ ' 



where 



= exp 



--(E/E,ay 



(16) 



(17) 



The value E^d is a function of E, B, p, T, jj. (see 
Ivan Adelsberg & Laf200 6' for details). 

The solution of the radiative transfer equation (|5]l is checked 
for the energy balance equation (|9]l, together with the surface 
flux condition 



47r 



0) dv = o-Hjj. = AnHi 



10, 



(18) 



Jo 



where the Eddington flux at any given depth m is defined as 
Hv(m)^—Y, dip j ^ tilUfi,^,m)dfi. (19) 



1=1 

The relative flux error 
Ho 



enim) = 1 



nOO ' 

Jq Hy(m)dv 



(20) 
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Fig. 5. Emergent spectra and temperature structures of fully ion- 
ized hydrogen model atmospheres of neutron stars with Teff - 
5 X 10^ K and log g - 14.3. The models with and without (dash- 
dotted curves) strong (B = 5 x 10'"^ G) magnetic field as well 
as with (solid curves) and without (dashed curves) vacuum po- 
larization effect are shown. The blackbody spectrum with the 
temperature equal the effective temperature is also shown in the 
upper panel (dotted curve). 

and the energy balance error as functions of depth 

-27T 
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dv 

Jo Jo 



(21) 



are calculated. 

Temperature corrections are then evaluated using three dif- 
ferent procedures. The first is the integral A-iteration method, 
modified for the two-mode radiation transfer, based on the en- 
ergy balance equation (|9|l. In this method the temperature cor- 
rection for a particular depth is found from 



AT. 



T'^vdv' 



where 



(22) 



(23) 



and a'y = o-y/(ky + crp. Averaged opacities are defined as 

1 pin r*+l 

= — j ^ kl(ti,<p)dij, (24) 

^vdiag diagonal matrix element of the A operator depend- 
ing on the mean optical depth in given mode The mean optical 

depth is defined by dr^ = {kl + (fpdm (see l Kuruc3ll97d for 
details of Avdiag - t"v dependence). This procedure is used in the 
upper atmospheric layers. The second procedure is the Avrett- 
Krook flux correction, which uses the relative flux error e//(m), 
and it is performed in the deep layers. And the third one is the 
surface correction, which is based on the emergent flux error (see 
Kurucz 1970 for a detailed description). 

The iteration procedure is repeated until the relative flux er- 
ror is smaller than 1% and the relative flux derivative error is 
smaller than 0.01%. As a result of these calculations, we obtain 
a self-consistent isolated NS model atmosphere, together with 
the emergent spectrum of radiation. 

Our method of calculation has been tested by a compari 



son to models for magnetized NS a tmospheres fSh ibanov et al 



199 1 iPavlov etaPl 1 994t IHo & Lail l200 1 : Ozel 2001l: lHo & Lai 
2003 t). Model atmospheres with partially ionized hydrogen are 
compared to models computed by Ho et al. (2007). We have 
found that our models are in a good agreement with these cal- 
culations. Our results are presented in Fig.|5j where the temper- 
ature structures and the emergent spectra of models with and 
without vacuum polarization are compared to a model with- 
out magnetic field. In Fig. |6] we present emergent spectra and 
temperature structures of the semi-infinite an d thin model atm o- 
spheres with the same parameters as used bv lHo et al] (l2007h . 

3. Results 

In this work we use the developed code mainly for studies of thin 
atmospheres above a condensed NS surface. In all calculations 
below we use the same surface gravity, \ogg - 14.3. 

One of the problem related to magnetars (in particular, 
AXPs) is the lack of any absorption feature at the proton cy- 
clotron energy, although early models of the magnetized NS 
atmo spher es predicted a strong feature at t he magnet a r field 
strengths (iHo & Laill20"oTt IZane et al1l200lh . IHo & Lail (l2003h 
suggested that a possible solution of this problem is the sup- 
pression of the cyclotron absorption feature due to the vacuum 
polarization. We confirm this result (see Fig.|5]l. In addition, we 
demonstrate that this absorption line is further reduced in a thin 
atmosphere without vacuum polarization. Figure|2]demonstrates 
emergent spectra of thin atmospheres without allowance for the 
vacuum polarization effect. In the semi-infinite atmosphere, a 
wide proton cyclotron line forms in agreement with the results 
of lHo&Lail (I200lh : IZane et al.1 (1200 lb . However, the absorp- 
tion feature disappears with decreasing the atmosphere surface 
density 2. The thin atmosphere is transparent to the continuum 
and absorption line wings, therefore the emergent spectrum ap- 
proaches the spectrum of the condensed surface everywhere ex- 
cept for a naiTow energy band at the center of the cyclotron line. 
The width of this absorption depends on the atmosphere thick- 
ness. 

So me XD INSs and CCOs sho w one or two absorption fea- 
tures (lHaberl 2007; Sanwal et al. 2002 ; Schwop e et al. 200^. 
Various hypothesis (so far inconclusive) were considered for 
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Fig. 6. 7b/? panel: Emergent spectra of the partially ionized hy- 
drogen semi-infinite (dashed curve) and thin (2 - l.lgcm^^, 
solid curve) model atmospheres (Teg - 5.25x10*' K, B - 4x10'^ 
G) together with the corresponding blackbody spectrum (dot- 
ted curve). The spectra are calculated with gravitational redshift 
(z=0.22) taking into consideration. Bottom panel: Temperature 
structures of the models from the top panel together with the 
temperature structure of a fully ionized hydrogen model (dotted 
curve). 

an explanation (see iMori & Hoi l2007l and references therein). 
Here we suggest another one, which we name "sandwich at- 
mosphere". A thin, chemically layered atmosphere above a con- 
densed NS surface can arise from accretion of interstellar gas 
with cosmic chemical composition. In this case, hydrogen and 
helium quickly se parate due to the high gravity (according to 
iBrown et alj|2002l the He/H stratification timescale can be es- 
timated as ~ p[-^ "'^^M ^' where pi = p/10 g cm~^, T(, - 
T/IQ^ K, and gi4 - g/lO'"* cm s"^). In this "sandwich atmo- 
sphere" a layer of hydrogen is located above a helium slab, and 
the emergent spectrum has two absorption features, correspond- 
ing to proton and ff-particle cyclotron energies. In Fig. [8] the 
emergent spectrum for one of such models is shown. The emis- 
sion feature at the helium absorption line arises due to a local 
temperature bump at the boundary between the helium and hy- 
drogen layers. This bump arises due to sharp changing of the 
plasma density and the opacity between helium and hydrogen 
layers. Clearly, some transition zone with mixed H/He chemical 
composition must exist between layers, and this rapid tempera- 



Fig. 7. Emergent spectra of the thin fully ionized hydrogen at- 
mospheres above a solid surface with T^ff = lO*" K and B - 10''* 
G in comparison to the semi-infinite atmosphere (dashed curve). 
The spectra of atmospheres with surface densities 'E-l (solid 
curve) and 100 g cm"^ (dash-dotted curve) together with the cor- 
responding blackbody (dotted curve) are shown. Vacuum polar- 
ization effect is not included. 



ture change can be reduced. We plan to calculate models with 
this kind of transition zone in future work. 

For comparison to observations, it is necessary to integrate 
the local model spectra over the NS surface. The effective tem- 
perature and magnetic field strengths are not uniform over the 
NS surface, and ge nerally the mag netic field is not perpendicular 
to the surface (see lHo et al.ll2008h . Therefore, it is necessary to 
compute model atmospheres with inclined magnetic field. This 
possibility is included in our code. For example. Fig. |9] shows 
spectra and temperature structures of model atmospheres with 
magnetic field perpendicular and parallel to the NS surface. 

Most of the XDINSs ha ve magnetic fi elds B > 10'^ G and 
color temperatures » 10^ K (Haberf2007'). Hydrogen model at- 
mospheres are partially ionized under these conditions and the 
vacuum polarization effect is also significant. Here we present 
first results of modeling of partially ionized hydrogen atmo- 
spheres using our radiative transfer code. In Fig. [TO] we com- 
pare spectra and temperature structures of the partially ionized 
hydrogen model atmospheres with and without the partial mode 
conversion effect. When the X-mode (having smaller opacity) 
partially converts to the O-mode in the surface layers of the at- 
mosphere, the energy absorbed by the O-mode heats these upper 
layers. As a result, the emergent spectra are closer to the black- 
body. 

For some of the XDIN Ss, optical counterparts have been 
found ( Mignani et al.ll2007l) . The observed optical fluxes are a 
few times larger compared to the blackbody extrapolation from 
X-rays t o the optical ran ge (see top panel of Fig. [12] for illus- 
tration). iHo et al.l (l2007h demonstrated, that a single partially 
ionized thin hydrogen atmosphere can explain this problem in 
the case of brightest isolated NS RXJ1856. 4-3754: the model 
fits well both, the observed optical flux and the X-ray spec- 
trum. RX J 185 6.4-3754 has very low puls ed fraction of radia- 
tion (a; 1.2 %. iTiengo & Mereghettill2007l) . therefore it is pos- 
sible to fit the radiation of this star by the single model atmo- 
sphere. Oth er XDINSs have larger pulse fractions, up to 18% 
(RBS 1223. lHaberl et ak 2004). In this case the temperature dis- 
tribution across the NS surface is not homogenous, and the ex- 
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Fig. 8. Emergent spectrum (top panel) and temperature struc- 
ture (bottom panel) of the "sandwich" model atmosphere above 
a solid surface with Teff = 10^ K and B - 10'^ G (solid curve) 
in comparison with the thin fully ionized hydrogen atmosphere 
(dashed curve) with the same parameters. The surface densities 
of both model atmospheres are 100 g cm in the "sandwich" 
model the H slab has 25 g cm"^ surface density and the He slab 
has 75 g cm"^. The corresponding blackbody spectrum (dotted 
curve) is also shown in the top panel. 



cess optical flux can be explained by the radiation from cool 
surface parts dSchwope et al.ll2005l) . 

We now investigate properties of partially ionized hydro- 
gen modes, which can be applied to the RBS 1223 atmosphere. 
The color temperature of this star, found from X-ray spectra fit- 
ting, is close to 10*^ K, with magnetic field B 4 x 10'^ G 
jSchwope et al.ll2007h . In particular, we investigate the optical 
flux excess in comparison to the X-ray fitted blackbody flux in 
this kind of models. For this aim we have calculated two sets of 
models with vacuum polarization and partial mode conversion. 
The models in the first set have effective temperatures T^ff -10^ 
K and the models of second one have effective temperatures 
Teff = 1.2 X lO*" K. In both sets B = 4 x 10'^ G, and models 
with surface densities 2 = 1, 3, 10, 30, 100 and 10^ (semi-infinite 
model) g cm"^ are computed. In Fig.[TT]we show emergent spec- 
tra and temperature structures for some models from the second 
set. Clearly, the X-ray spectra of the models with Z < 10 g cm - 
are close to a blackbody and, therefore, better fit the observed 
X-ray spectrum. 



> 




0.1 1 10 
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Fig. 9. Emergent spectra and temperature structures of partially 
ionized hydrogen model atmospheres with Teif - 1.5 x 10* K 
with inclinations of the magnetic field {B = 4x 10'^ G) to the sur- 
face normal 6b equal 0° (solid curves) and 90° (dashed curves). 
The vacuum polarization eff'ect is not included. The correspond- 
ing blackbody spectrum is also shown in the upper panel (dotted 
curve). 

In Fig. [12] (bottom panel) we show the ratio of the model at- 
mosphere flux to the X-ray fitted blackbody flux in the optical 
band depending on Z for both sets. The observed ratio is about 
5 ( Kaplan et al. 2002b), in agreement with the semi-infinite at- 
mosphere models. However, the observed blackbody like X-ray 
spectrum agrees with the thin atmosphere models, for which 
this ratio is close to 1. Therefore, the observed optical excess 
cannot be explained by the thin atmosphere model alone; in- 
stead, it can arise due to a nonuniform surface temperature dis- 
tribution, in agreemen t with the RBS 1223 light curve modeling 
(ISchwope et al.ll2005l) . 

4. Conclusions 

In this paper a new code for the computation of magnetized NS 
model atmosphere is presented. It can model fully ionized and 
partially ionized hydrogen atmospheres in a wide range of effec- 
tive temperatures (3 x 10^ - 10' K) and magnetic fields (10'^ - 
10'^ G), with any inclination of the magnetic field to the stellar 
surface. The vacuum polarization effect with partial mode con- 
version is taken into consideration. Calculated emergent spectra 
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Fig. 10. Emergent spectra and temperature structures of the par- 
tially ionized hydrogen model atmospheres with T^ff = 1.5 x 10'' 
K with (solid curves) and without (dashed curves) vacuum polar- 
ization effect (with the partial mode conversion). The magnetic 
field strength is B = 4 x 10'^ G. The corresponding blackbody 
spectrum is also shown in the upper panel (dotted curve). 



and temperature structures of the model atmospheres agree with 
previously published ones. 

We presented new results obtained using this code. We have 
studied the properties of thin atmospheres above condensed NS 
surfaces. We demonstrated that the proton cyclotron absorption 
line disappears in the thin hydrogen model atmospheres. A new 
thin "sandwich" model atmosphere (hydrogen layer above he- 
lium layer) is proposed to explain the occurrence of two absorp- 
tion features in the observed X-ray spectra of some isolated NSs. 

We analyzed the optical excess (relative to the X-ray fitted 
blackbody flux) in the model spectra of partially ionized hy- 
drogen atmospheres with vacuum polarization and partial mode 
conversion. 

A set of model atmospheres with parameters (effective tem- 
perature and the magnetic field strength) close to the probable 
parameters of the isolated NS RBS 1223 were calculated. We 
found the optical flux excess x 5 for the semi-infinite model 
atmospheres decreases down to 1 with decreasing surface den- 
sity E of the atmosphere. Spectra of thin model atmospheres fit 
the observed RBS 1223 X-ray spectrum better, therefore we con- 
clude that the observed optical excess should be explained by 
nonuniform surface temperature distribution. 



Fig. 11. Emergent spectra and temperature structures of the par- 
tially ionized hydrogen model atmospheres with Tgff - 1 .2 x 10^ 
K with vacuum polarization effect the partial mode conversion 
shown for various surface densities E (solid curves - 1 g cm"^, 
dashed curves - 10 g cm"^, dash-dotted curves - semi-infinite at- 
mosphere). The magnetic field strength is B = 4 x 10'^ G. The 
corresponding blackbody spectrum is also shown in the upper 
panel (dotted curve). 

The accuracy of the thin and sandwich model atmospheres 
is currently limited by the inner boundary condition for the ra- 
diation transfer equation. We used blackbody radiation as this 
condition, but a higher accuracy can be achieved by replacing i t 
by the condensed surface condition dvan Adelsberg et ani2005h . 
This will be done in a future work. We are also planning to in- 
clude the effect of magnetic field and temperature distributions 
over the stellar surface to compute an integral emergent spec- 
trum from isolated NSs. 
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